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COMPUTING  THE  GRAZING  ANGLE  OF  SPECULAR  REFLECTION 


INTRODUCTION 

'  / 

In  this  report  we  give  three  methods,  for  computing  the  grazing  angle  of  specular  reflection.  The 
first  was  developed  by  Fishback  [1,2]  in  1943  and  refined  somewhat  by  Blake  [3]  in  1980,  and  jb  gives 
an  approximation  to  the  grazing  angle.  The  second  is  an  iterative  method,  whose  derivation  is  easy  and 
provides  for  the  computation  of  the  grazing  angle  to  any  degree  of  accuracy  specified.  The  third 
method  shows  that  there  is  a  closed  form  explicit  expression  for  the  grazing  angle.  This  in  turn  also 
allows  for  a  computation  to  an  arbitrary  degree  of  accuracy.  A  FORTRAN  program  of  the  second 
method  is  includedin  Appendix  B. 

GENERAL  REMARKS 

The  geometry  of  the  spherical-earth  problem  and  definitions  of  various  angles  and  distances  are 
given  in  Fig.  1. 

The  problem  of  computing  the  grazing  angle  of  specular  reflection  is  that  of  computing  $  in  Fig. 
1,  i.e.,  finding  ^  an  angle  between  the  tangent  line,  f,  to  the  circle  and  AR  where  <fi  is  also  the  angle 
between  land  TR. 
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We  assume  that  <p  -  (<A|  +  02);  h\,  A2,  and  rt  are  given.  Define  k, 


METHOD  ONE 


In  method  one.  let 


r,<P. 


r,(/ti  +  A2)  +  -7s2 

4 


and 

q  -  sin_l|2p_3r(s(/i2  -  A|)|. 

Then  when  h{  and  h2  are  very  much  less  than  r,,  Fishback  approximates  <f>\  by 

.  .  1  ,  p 

<P  1  -  —  sm 

*  rt 

Finally,  the  grazing  angle  ip  is  given  by 

ip  -  tan_1(cor  -  A,  esc  ^j). 


(1) 

(2) 

(3) 

(4) 

(5) 


The  error  in  the  approximation  for  <£,  (and  consequently  for  ip)  is  not  generally  known.  How¬ 
ever,  if  hi  =  A2,  then  <£1  ”  1/2 0  and  Eq.  (5)  provides  the  exact  result  for  the  grazing  angle. 

METHOD  TWO 

Method  two  uses  an  iterative  procedure  for  computing  the  grazing  angle  to  any  degree  of  accu¬ 
racy 


First,  applying  the  law  of  sines  to  triangles  OAR  and  ORT,  we  obtain 
</>i  +  <l>  -  cos-1  ( kx  cos  ip) 


and 


<t>2  +  “  COS  1  (,k2  cos  Ip). 


(6) 

(7) 


Adding  the  equations  (and  recalling  that  <f>  -  <f>i  +  <f>2)  we  have. 


Ip  -  g(</<) 

(8) 

where. 

gty)  -  y|cOS“'(A|  cos  «//)  +  cos-1  (Ar2  cos  Ip)  -  0J. 

(9) 

If  we  choose  0n  arbitrarily  and  define 

<Pi+ 1  ”  «(0,).  i  “  0,  1,  2,  ... 

(10) 

then  by  the  results  of  Appendix  A, 

lim  1 pi 

(11) 

exists  and  is  the  unique  solution  of  Eq.  (8),  i.e.,  (11)  is  the  grazing  angle. 
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Approximations  to  the  grazing  angle  are  given  by  the  terms  of  the  sequence  i/>0-  0i . with  suc¬ 

cessive  terms  providing  more  accurate  approximations.  In  practice,  when  the  relative  difference  of  suc¬ 
cessive  terms  of  this  sequence  differ  in  absolute  value  by  no  more  than  a  predetermined  constant,  we 
obtain  the  grazing  angle  to  our  required  degree  of  accuracy.  A  FORTRAN  program  implementing 
method  two  is  given  in  Appendix  B. 

To  obtain  a  rapid  convergence  for  method  two,  we  might  use  method  one  to  obtain  the  initial 
value,  <ji0-  In  addition,  method  two  may  be  useful  in  real-time  computation  of  the  grazing  angle,  since 
an  angle  once  computed  may  be  used  as  the  initial  value  for  an  update  computation. 

METHOD  THREE 

The  third  method  produces,  in  principle,  an  explicit  expression  for  the  grazing  angle. 

Let  U  -  exp(/<£)  and  Z  -  exp(2/«/») - 

Then  replacing  U  and  Z  in  Eq.  (8),  we  obtain  the  following  quartic  (see  Appendix  C  for  deriva¬ 
tion): 

aZ4  +  /SZ3  +  CZ2  +  0Z  +5-0,  (12) 

where, 

A  -  U-  kxk2, 
a  —  UA, 

0  —  k\  +  k2  —  2kxk2U, 
and 

C-2Re[/3-  Ua[. 

Since  Eq.  (12)  is  a  quartic,  the  roots  can  be  exhibited  explicitly  using  the  classical  method  of  Fer¬ 
rari  and  Cardan  [4],  At  least  one  of  these  roots  lies  on  the  unit  circle.  Let  Z.  designate  any  of  those 
roots  on  the  unit  circle.  Then 

s  y  cos"1  {Re  (Z.)},  (13) 

and  that  unique  value  of  t/>.  that  satisfies  Eq.  (8)  is  the  grazing  angle. 

Rather  than  use  the  method  of  Ferrari  and  Cardan  to  find  the  four  roots  of  Eq.  (12),  it  is  easier 
and  more  efficient  *o  solve  the  quartic  numerically  on  a  computer  using  a  polynomial  root  finder  rou¬ 
tine. 

CONCLUSION 

Three  methods  for  computing  the  grazing  angle  of  specular  reflection  are  given.  The  first  pro¬ 
vides  an  approximation  where  the  error  is  not  known.  Methods  two  and  three  will  provide  computa¬ 
tions  good  to  any  degree  of  accuracy.  Method  two,  an  iterative  procedure,  may  be  especially  useful  in 
real  time  computation. 
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Appendix  A 

PROOF  OF  CONVERGENCE  TO  THE  GRAZING  ANGLE 


Let  X  be  an  arbitrary  set,  d  a  metric  on  X,  and  /  a  contraction  mapping  of  X  to  itself,  if  is  called 
a  contraction  mapping  if  there  is  a  positive  constant  k  <  1  such  that 

difix),  fiy))  <  k  dix,y)  for  all  x,  y  in  X). 

The  following  result  may  be  found  in  [All. 

Theorem:  Every  contraction  mapping,  /,  on  a  complete  metric  space  has  a  unique  fixed  point, 
i.e.,  there  is  a  unique  x  in  X  such  that  fix)  -  x  Furthermore,  for  an  arbitrary  x0  in  X,  the  sequence 
given  by  x„+1  -  fix,),  n  -  0,  1,2...  converges  to  x 

As  a  special  case,  let  X  be  the  reals  and  for  x,y  in  X  let  dix.y)  -  |x  -  y  I.  Then  iX,d)  is  a  com¬ 
plete  metric  space. 

Now  let  /  be  a  differentiable  function  defined  on  X  and  suppose  furthermore  that  there  is  a  con¬ 
stant  k,  0  <  k  <  1  such  that  |/*(x)|  <  k  for  all  xin  X.  Then  /is  a  contraction  mapping.  This  may  be 
seen  as  follows.  Let  x  and  y  be  real  numbers,  x  <  y.  By  the  mean  value  theorem  there  is  a  number  (, 
x  <  f  <  y  such  that 

fix)  —  fiy)  **  (x  —  y)f(f)' 

Hence, 

I  fix)  -  fiy)  I  -  lx  -  I !/’(()  I  <  *  lx  -  y\. 

These  results  may  be  summarized  in  the  following. 

Corollary:  Let  A  be  the  reals,  k  a  constant  0  <  k  <  1,  and  /  a  differentiable  function  on  R 
such  that  l/*(x)  I  <  k  for  all  x  in  R.  Then  there  is  a  unique  real  number  a  such  that 

a  -  fia). 

Moreover,  if  x0  is  an  arbitrary  real  number  and  x.+t  -  fix,),  n  -  0,  1,  2 . then 

lim  x,  -  a. 


With  respect  to  Method  2,  we  have  given  the  differentiable  function 

gix)  -  y[cos-,(*i  cos  x)  +  cos-,(*2  cos  x)  -  ^j. 

Differentiating,  we  have 

k,  sin  x  k2  sin  x 

2  g'ix)  -  ....  j  +  -7=  ,  ...  , 

Vl  -  kf  cos 1  x  Vl  ~  k]  cos?x 
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Recalling  that 


k,  - 


rf  +  h, 


<  1, 


1,  2 


we  have 


|2g'(x)l  ^  k 


Vi  -  cos2  x  r 

1  -  k\  cos2  X  2V  1 


1  —  COS2  X 


—  k]  cos2  x 


4  k,  +  k2  <  2. 

or 

l*'(x)|  <  ~~2—  <  1. 

Hence,  by  the  Corollary,  there  is  a  unique  number  a  such  that 

a  -  g(a). 

Moreover,  if  x0  is  arbitrary  and  xB+1  -  g(xn),  n  -  0,  1,2,...,  then 

lim  x„  —  a. 
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Appendix  B 

A  FORTRAN  PROGRAM  OF  METHOD  TWO 


SOURCE  LISTING  ASC  FAST  FORTRAN  CONPILER  RELEASE  FTF»05?9.P29AIRC 

ST4TE-ENT  CP  OPTIONS  «  <•*,*)  OATE  =  05 /06/82T  8  2 . 1  2f>  ) 


SUBROUTINE  GRATEf Hl.H?, THETA.flE.PHI> 

r 

C  THIS  POUT  INF  CONPUTES  THE  GRATING  ANGLE  PHI  POP  SPHFRtf.AL 

C  M»TH  SPECULAR  REFLECTION. 

C 

C  m-=  ANTENNA  HEIGHT. 

C  H2=  TARGET  HEIGHT. 

C  OF-  EARTH  EFFECTIVE  RADIUS  OR  EARTH  RADIUS* 

C  THETA  IS  CENTRAL  ANGLE  IN  RADIANS. 

c 

C  PHI  «U$T  BE  INITIAL  I TfO  IN  CALLING  ROUTINE:  PHI  .NF.  0. 

C  PHI*  GRATING  ANGLE  OUTPUT  IN  R  AO  I AN  $ . 

C  RELATIVE  ERROR  IN  phi  TS  10*«-R.  RUT  CAN  BE  DECREASED  B»  CHANGING 

C  VALUE  OF  TflL  IN  DATA  STATE  RENT  . 

C 

C  TNPUTS  ANO  OUTPUTS  ARE  IN  0£ AL* 8 1  RE AL*A  CAN  RE  USED  OY  RENOVING 

C  ILLICIT  STATENENT  AND  ADJUSTING  TflL. 

r 

C  Hi,  HJ,  RE.  THETA  are  inputs.  UNITS  for  Mj,  M2,  re  «ust  be 
C  CONSISTENT. 

C 

c 

INPLTCIT  RFAL*8<A-M,0-I) 

c 

DATA  TOL/l.O-8/ 

C 

RKI«RE/<HIPRE) 

»X2*RE/<N2*RE) 

C 

A=OARCOSCR*I>»OARCOS<RR2>-THFTA 

JF(  A  .GE.  0.00  .AND.  A  .LE.  1.0-lS  )  GOTO  U 

C 

20  GsO.ROOPfOARCOSCRKlPnCOSTPHI )  T»D ARCOSCRK2 *0COS<PHT > I-THET A ) 

C 

RTST*CG-PHI)/PHI 

®TST*0*8S<RTST> 

C 

IFCRTST.LE.  TOL  >  GOTO  10 

r. 

PMT=r, 

c 

GOTO  20 

c 

10  PMT=G 
C 

RETURN 

f. 

It  PHT»0.O0 
RETURN 
f  NO 
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Appendix  C 

DERIVATION  OF  THE  QUARTIC 


Beginning  with  Eq.  (8),  we  have 

cos  (2 0  +</>)“  k]k2  cos2 1 h  —  —  k 2  cos2 1 p  1  —  k2  cos2 1 It  . 

Squaring,  we  obtain 

(1  -  A:2  cos2  0)0  -  k2  cos2 1 /,)  =  [Art A2  cos2 1 1/  -  cos(2i p  +  0)]2, 
or 

1  —  (k2  4-  k2)  cos2 1 li  —  cos2(2i ^  +  0)  —  2A| A2  cos2  0  cos(20  +  0).  (Cl) 

U  —  exp  (/0)  and  2  —  exp  020), 


cos2  0  ■*  -j-(2  +  2  +  2), 

4 


cos  (20  +  0)  -  j[i/2  +  UZ],  (C3) 

Multiplying  both  sides  of  Eq.  (Cl)  by  4 U2Z2  and  substituting  Eq.  (C2)  and  Eq.  (C3)  into  Eq. 
(Cl)  we  have 

4 U2Z2  -  U2(k?  +  k})(Zl  +  2 Z2  +  2) 

-  (U*Z*  +  2 U2Z2  +  1)  -  Uk\k2(Z2  +  2 Z  +  1  )W2Z2  +  1). 

Then  grouping  the  terms  in  powers  of  2  gives 

U2AZ*  +  l/2/S23  +  £/2( 2Re  C/3  -  t/I ))  Z2  +  U2pZ  +  UA  -  0. 

Now  dividing  by  U 2  gives  Eq.  (12). 


Let 

so  that 

and 
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